Dispersion control for matter waves and gap solitons in optical superlattices 
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We present a numerical study of dispersion manipulation and formation of matter-wave gap 
solitons in a Bose-Einstein condensate trapped in an optical superlattice. We demonstrate a method 
for controlled generation of matter-wave gap solitons in a stationary lattice by using an interference 
pattern of two condensate wavepackets, which mimics the structure of the gap soliton near the edge 
of a spectral band. The efficiency of this method is compared with that of gap soliton generation in 
a moving lattice recently demonstrated experimentally by Eiermann et al. [Phys. Rev. Lett., 92, 
230401 (2004)]. We show that, by changing the relative depths of the superlattice wells, one can 
fine-tune the effective dispersion of the matter waves at the edges of the mini-gaps of the superlattice 
Bloch-wave spectrum and therefore effectively control both the peak density and the spatial width 
of the emerging gap solitons. 
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I. INTRODUCTION 

In recent years Bose-Einstein condensates (BECs) 
loaded into optical lattices have become an important 
tool in the studies of linear and nonlinear behavior of 
coherent waves in periodic systems. The band-gap struc- 
ture of the matter-wave spectrum resulting from the pe- 
riodicity of the optical trapping potential is responsible 
for much of the observed coherent behavior of the con- 
densate, including Bloch oscillations and Landau-Zener 
tunnelling 1', ^ . One of the more dramatic effects of the 
lattice on the BEG dynamics is the dependence of the 
group-velocity dispersion of matter waves on the curva- 
ture of the spectral band. The latter, in turn, is a func- 
tion of the wavepacket quasi-momentum in the lattice 
rest frame. Due to the possibility to create a moving 
optical lattice by varying the relative detuning of the in- 
terfering laser beams, the condensate can be loaded into 
the lattice with a well-defined quasi-momentum . This 
enables precise engineering of the dispersion of the BEG 
wavepackets in moving optical lattices 0, Q . 

Dispersion management of a BEG wavepacket loaded 
into a one-dimensional (ID) moving lattice was recently 
investigated experimentally ^ |M 1^ > in both linear (non- 
interacting) and non- linear (weakly interacting) regimes. 
In the nonlinear regime, a balance between the effects of 
nonlinearity and dispersion can produce a bright matter- 
wave soliton, a localized BEG wavepacket that maintains 
a constant spatial structure. The creation of a bright 
soliton in a repulsive condensate (as opposed to an at- 
tractive condensate) is only possible due to the negative 
effective dispersion at the edge of the first Brillouin zone 
(BZ) of a ID optical lattice (see e.g. Q)- Driving the 
condensate from the middle to the edge of the Brillouin 
zone achieves transition between the regimes of positive 
and negative effective dispersion. Evolution of the matter 
wavepacket in the negative dispersion regime can result 
in the self-focusing of the repulsive BEG in the form of a 
fundamental gap soliton 0. 

The absolute value of the dispersion experienced by 



a BEG wavepacket can be varied by changing the lat- 
tice depth [3 AH previous experimental studies on 
matter- wave dispersion management and gap solitons in- 
volved single-periodic shallow lattices P, ^ Q Q which 
are characterized by narrow spectral gaps, greater cur- 
vatures of the spectral bands, and hence larger values 
of the effective dispersion at the gap edges. Stronger 
dispersion requires larger matter-wave nonlinearity, and 
hence larger atomic density, and/or wider wavepackets 
to achieve the localized state Therefore, dispersion 
manipulation in an optical lattice can potentially deliver 
control over the characteristics of BEG solitons. 

In this paper, we explore novel possibilities to con- 
trol the magnitude of dispersion experienced by a BEG 
wavepacket at the edges of the spectral hands by modi- 
fying the shape of a double-periodic optical superlattice. 
As known from the theory of diatomic lattices and semi- 
conductor superlattices, an extra periodicity opens up 
additional narrow stop-gaps (or mini-gaps) in the band- 
gap spectrum. Here we show that the values of effective 
dispersion at the edges of these mini-gaps can be var- 
ied within a much greater range than it is possible for a 
single-period lattice of a reasonable depth. 

To demonstrate the efhcient dispersion control and 
generation of immobile gap solitons, we investigate the 
process of nonlinear localization in both moving and sta- 
tionary superlattices. In the former case, in order to ac- 
cess the gap regions, the initially stationary condensate 
in the middle of the ground-state Brillouin zone is driven 
to the band edge by accelerating the lattice. By the 
end of the adiabatic acceleration process, the condensate 
wavepacket has a quasi-momentum and an internal struc- 
ture of a linear Bloch wave at the gap edge. In this paper, 
we numerically simulate this process, starting from the 
initial adiabatic loading of the condensate into a quasi- 
ID superlattice. Our numerical simulations closely follow 
the current experimental procedure for matter-wave gap- 
soliton generation 

In addition, we suggest a new method for creating the 
correct initial conditions for the soliton formation at the 
gap edge in a stationary lattice, by interfering two identi- 



cal wavepackets with equal momenta of the opposite sign, 
corresponding to the opposite edges of the first BZ. This 
enables us to efficiently create a matter wavepacket with 
the correct internal structure at the relevant gap edge 
and avoid large time scales associated with the adiabatic 
acceleration. Such a wavepacket preparation technique 
was previously explored in optics Here 
we show that this method leads to more efficient soliton 
generation and shape control of the emerging gap solitons 
via dispersion control in the superlattice. 



II. MATTER- WAVE SPECTRUM IN A 
SUPERLATTICE 



In the mean-field approximation, the dynamics of 
a Bose-Einstein condensate is described by the Gross- 
Pitaevskii (GP) equation. Assuming that the condensate 
cloud is elongated in the direction x, with the ratio of 
the corresponding frequencies at most = w^/wj^ ^ O.I 
[M IM> Q| I the condensate wavefunction in the axial di- 
mension can be described b y t he one-dimensional (ID) 
GP equation (see, e.g., Ref. jla|): 
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where V{x^ t) is the external trapping potential, and giD 
characterizes the strength of the two-body interactions 
rescaled for the case of the ID geometry. The structure 
of the condensate wavefunction in the transverse dimen- 
sions is determined by a tight harmonic potential. 

Assuming that any additional trapping along the axial 
direction (e.g., due to a magnetic trap) is weak, we write 
V{x,t) as a periodic potential of a ID optical superlat- 
tice: 

V{x) = U[esm^{Kix) + [l - e)iiin^{K2x)]. (2) 

The superlattice potential given by Eq. (|2Jl can be ob- 
tained by superimposing two independent (either de- 
tuned from each other T^l or orthogonally polarized 
far off-resonance single-periodic standing waves with dif- 
ferent periods di and ^2- The larger of the two periods, 
e.g. d = di, defines the coarse periodicity of the lattice. 
The lattice wavevectors are given by Ki = n/di and 
K2 = 7r/d2 where the commensurable periods are chosen 
such that di/d2 = 2. The relative and total intensity 
of the standing waves are controlled by the parameters 
< e < I and U respectively. 

Equation with the lattice potential given by 

(PI, was made dimensionless by using the characteris- 
tic length fli — d/iT, energy El — h? /ma\ = 2E,-cc, 
and time cj^ — H/El, where m is the atomic mass. 
In these dimensionless units, and with the original 
three-dimensional condensate wavefunction normalized 
by a^'^^^, the nonlinear coefficient in Eq. (I) becomes: 
giD = 2(as/aL)(wj^/wL). Here we use the parameters 
of the ^''Rb condensate: m = 1.44 x I0~^^ kg and 




FIG. 1: The structure of the superlattice potential described 
by Eq. (|5J for different values of e and a fixed potential height 
Vb = 1: single-periodic optical lattice for (a) e = and (d) 
e = 1, and double-periodic superlattices for (b) e = 0.3 and 
(c) e = 0.7. 



as — 5.7 nm. We assume that the periods of the standing 
waves forming the superlattice are d = di = 350 nm and 
d2 — 175 nm, so in the dimensionless units, Ki ~ I 
and K2 = 2. With stron g tr ansverse confinement of 
cj_L 500 - 550Hz H El 111, giD ~ 0.001, and this 
is the value we use for most of the simulations in this 
paper. 

The shape of the optical superlattice Q depends on 
the values of the parameters, U and < e < I. In the 
limits e ^ or £ ^ I, the lattice becomes single-periodic, 
and U coincides with the height of the lattice Vq. For 
e 7^ and e ^ I, Eq. Q describes a double-periodic 
superlattice with U defined through the amplitude of the 
periodic potential Vq as 



U = I6Fn 



(1-g) 
(4-3e)2- 



Figure ^ shows several examples of the superlattice po- 
tential (01 with the constant amplitude Vq — 1. As seen 
in the figure, the relative depth of the large and small lat- 
tice wells can be manipulated by varying e, whilst keeping 
the height of the lattice constant. 

In the linear regime {giD = 0), the stationary solutions 
to Eq. with the periodic potential are Bloch waves: 



il}{x, t) = (j){x)e 



u{x)e 



c , 



(3) 



where ^ is the matter-wave chemical potential and the 
function u{x) has the periodicity of the lattice. The lin- 
ear matter-wave spectrum /i(fc) consists of bands where 
the real wavenumbers fc(/i) correspond to the oscillatory 
Bloch- wave solutions. The bands are separated by "gaps" 
where Im(A:) 7^ 0. Figure |2Ia) shows the band-gap dia- 
grams plotted in the extended zone scheme for a single- 
periodic lattice (e = 0) and a superlattice (e = 0.05), 
both for the lattice amplitude of Vb = 5. Since the coarse 
periodicity of the superlattice is twice that of the single- 
periodic lattice, the size of the Brillouin zones of the su- 
perlattice are half that of the single-periodic lattice at 
e = 0. This leads to mini-gaps appearing in the super- 
lattice band-gap spectra at A: = I. A shaded stripe in 
Fig. |2Ia) shows the lowest energy mini-gap at e = 0.05. 
As e grows, the size of the mini-gap increases. 
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FIG. 2: (a) The band-gap spectra for a single-periodic optical 
lattice (e = 0) and a double-periodic superlattice with weak 
modulation (e — 0.05). In both cases, the lattice height is set 
to Vb = 5. The first mini-gap is shaded, (b) The spectrum 
around the first mini-gap for a strongly modulated superlat- 
tice at Vb = 1 and e = 0.3. 
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FIG. 3: The effective group-velocity dispersion D (a) at the 
edges of the first gap of a single-periodic optical lattice of 
varying strength, and (b) at the edges of the first mini-gap 
for a superlattice of height Vb = 1 and period d = tt, as the 
shape of the superlattice is changed by varying e. 

For an interacting condensate {giD =/= 0), the solu- 
tions of the model equation near the band (or mini- 
band) edges can be sought in the form: tp{x,t) = 
exp {—iijLt)f{x,t)(f){x), where (f>{x) is the linear Bloch 
wave at the correspondin g ba nd edge and f{x,t) is a 
slowly varying envelope I3,[i3)l23- Then, the dynamics 
of the envelope are governed by the reduced GP equa- 
tion 0, El 113: 

df(x,t) ( D d"^ . , ,,2! , 

where D — d^^/dk^ = dvg/dk is the effec- 
tive group velocity (vg) dispersion and gm — 
giD / |0(a;)|'*dx/ J \(j){x)\'^dx. 

Figure |3Ia) shows the effective group velocity disper- 
sion, D, near the lowest energy gap of a single-period lat- 
tice for varying lattice height. As seen in Fig. (SJa), the 
dispersion is negative at the bottom of the gap and posi- 
tive at the top. As the condition gioD < is required for 
the formation of bright gap solitons, bright gap soliton 
families originate near the bottom of the gaps, for repul- 



sive condensates {gio > 0), and near the top of t he gaps , 
for attractive condensates {gio < 0) pl ll2a. l2l. I m l25l| . 
and exist for the entire range of chemical potentials 
within that gap |2^ . The group velocity Vg — 

dfi/dk vanishes at the band edges, hence the gap soli- 
tons form as immobile localized wavepackets. We note 
that formation and dynamics of gap solitons could be 
significantly modified in the case when the harmonic con- 
finement in the longitudinal dimension becomes impor- 
tant I23, H^l- However, provided the longitudinal 
harmonic trapping is weak, the band structure imposed 
by the lattice is still well defined. Moreover, in many 
experiments the longitudinal trapping potential is either 
sufficiently flat on the scale of hundreds of lattice wells, 
or is removed altogether after the condensate is loaded 
into an optical lattice 

Gap soliton formation can also occur in the minigaps 
of the superlattice spectrum. The size of the mini-gap 
increases with growing e. This decreases the curvature 
in the band structure immediately surrounding the mini- 
gap, leading to a decrease in the magnitude of the effec- 
tive dispersion at the gap edges. Figure |2Ib) shows that 
simple variation of the relative depths of the two wells 
in the superlattice (i.e. e) at a constant lattice height 
(Vb = 1) opens up access to a large range of effective 
dispersions at the edges of the mini-gap. As Figure |2|[a) 
shows, this kind of dispersion control can also be per- 
formed in a single-periodic lattice by changing the lat- 
tice height. However, the range of effective dispersion 
values in a single-periodic lattice is much smaller than 
that achievable in a superlattice, and the larger propor- 
tion of this range is only accessible for Vb < 1. 



III. GAP-SOLITON GENERATION IN A 
MOVING LATTICE 

The key to producing a spatial gap soliton in a BEG 
is to access the band-gap edge of the linear spectrum, 
where the suitable effective dispersion can assist spatial 
self-focusing of the coherent matter waves. For a repul- 
sive condensate, this occurs at the bottom edge of the 
first gap, where the dispersion is negative. In experi- 
ments, the stationary condensate is initially loaded into 
the lattice ground state by slowly ramping up the lat- 
tice height. This produces a wavepacket in the lowest 
energy band (n = 0) in the middle of the first Bril- 
louin zone (centered at fc = 0). To ensure that other 
energy bands are not populated, the ramping process 
should be adiabatic. To adiabatically load the conden- 
sate into the lattice at quasi- momentum fc, the change 
of the single-particle BEG Hamiltonian on the time scale 
of a single Rabi oscillation associated with the transition 
between the Bloch states in the ground state, 4>o,k, and 
the first excited band, (/)i,fc, should be much smaller than 
the energy gap between the n = and n — \ bands 
at quasi-momentum k. In our dimensionless units, this 
adiabaticity condition translates to AVo/^r ^ A/i^(fc) 



4 



1600 

t 

1200 


(a) 


''' , 


800 






400 




I 


n 







(b) 

1 


f=320 

■J 




(c) 

ill 


1 



-400 -300 -200 
X 

1600 



-100 -450 -300 -150 




(e) 


« = 




..illNili.. 



(f) 


t = 1600 


.,11 


illillilL.. 



-200 -100 



100 200 -150 



FIG. 4: Soliton generation by (a) a single Gaussian 
wavepacket accelerated to the edge of the first mini-gap 
{k = 1) in a moving superlattice and (d) interference of two 
counter-propagating Gaussian wavepackets in a stationary su- 
perlattice. Both superlattices have s = 0.3 and Vo — 1. Fig- 
ure (a) is plotted in the lattice frame moving at iircc. Initially 
the wavepacket contains 500 atoms, and has a width param- 
eter of w = 100. In (d) the initial parameters are w = 81.6 
and N = 500. (b,c) and (e,f) capture the spatial structure of 
the wavepackets at different evolution times. 



(S m 113), where is the duration of the lattice 
ramp-up, AVq is the total change of the lattice ampli- 
tude, and A/i(fc) is the size of the first (mini)gap at 
the quasi-momentum k. In the middle of the Brillouin 
zone (BZ), a,t k — 0, the size of the first mini-gap is 
equal to four single-photon recoil energies 0. There- 
fore, the single-particle adiabaticity condition for load- 
ing a stationary condensate into the lattice at fc = is, 
AI/oAr< A/i2(fc = 0) = 4. 

In order to access the band edge, a constant accel- 
eration is introduced suddenly by applying time- 
dependent frequency detuning to the lasers used to create 
the lattice. The superlattice potential given by Eq.® be- 
comes time-dependent: V{x,t) — U[esm^{Ki{x — St)) + 
(1 — e)sm^{K2{x — 6t))]. Then, the reference frame in 
which the lattice is stationary (the lattice rest frame) is 
moving in the laboratory frame with velocity v{t) = S{t). 
If the acceleration, a = dS/dt, is constant, a constant and 
uniform force is applied to the condensate which is 
somewhat analogous to the application of a constant ho- 
mogeneous electric field to electrons in a crystal. The 
quasi-momentum of the condensate depends on the du- 
ration of the acceleration, ^a, as k{tA) = k{0) + atA, 
and in order to move the condensate to the band-edge 
{k{tA) = 1), the recoil velocity v^-cc = 1 of the lattice is 
reached with acceleration a = l/tx- The acceleration is 



adiabatic, i.e. the upper bands are not populated, pro- 
vided the probability of transitions between the n — Q 
and n — 1 bands are low. This requires a small mag- 
nitude of the acceleration, |A(5/iA| ^ A/i(fc) For 
acceleration at a constant lattice height, the critical en- 
ergy gap is at the edge of the Brillouin zone {k = ±1), 
where the gap is generally much smaller than that at 
A: = 0, as can be seen in Fig. |21 The ramping adiabatic- 
ity condition is also much harder to fulfil at the band-gap 
edge than in the middle of the Brillouin zone. This is why 
it is preferable to load the condensate into the lattice at 
A: = and then accelerate it to the band-edge instead of 
simply ramping up an already moving lattice at velocity 

We model the process of ramping and accelerating the 
lattice by starting with an initially stationary Gaussian 
wavepacket: il){x,t = 0) = ^exp(— x^/w^) where A is 
the amplitude and w is the width. Figure EJa) shows 
the three stages of the evolution of an initially Gaussian 
wavepacket with w = 100 and TV — 500 atoms in a super- 
lattice with e = 0.3 and a potential depth of Vq = 1 [the 
band-gap structure of this lattice is shown in Fig.|2Ib)]: 

• The superlattice is ramped up from Vq = to 
Vo 1 in time = 20. This places the wavepacket 
at fc = in the lowest energy band [point A in 
Fig. I2Ib)]. The wavepacket develops an internal 
structure with the same spatial properties as a 
Bloch wave at A: = 0. 

• The superlattice is accelerated to the edge of the 
first mini-gap [point B in Fig. |2Ib)] in time tp^ = 
300. 

• The wavepacket is then left to evolve in a super- 
lattice moving with constant velocity Wrcc i-e- at 
the edge of the Brillouin zone. Figures Ela-c) are 
plotted in a frame moving at the recoil velocity Urcc 
so that the wavepacket appears stationary when it 
is at the gap edge where its group velocity with 
respect to the lattice is zero. 

For sufficiently long times of the wavepacket evolution 
at the mini-band edge, a localized wavepacket with the 
spatial structure of a gap soliton emerges, as seen in 
Fig. E[c). A fundamental gap soliton near a gap edge 
takes the form of a broad sech-shaped envelope with 
an internal structure resembling the linear Bloch waves 
at that particular edge. This spatial structure allows 
sufficiently broad gap solitons to be described using a 
wave-envelope approximation with an appropriately cho- 
sen carrier wave, as described in the previous section. 
Due to large differences between the ideal shape of the 
soliton and the shape of the initial wavepacket, signifi- 
cant oscillations in the shape of the wavepacket arc ob- 
served [see Fig. 0Ja)]. This is also a well known part of 
soliton formation in the absence of the lattice [sil Is^ . 
The closer the initial width of the wavepacket is to that 
of the soliton, the more efficient is the process of soli- 
ton generation. This is illustrated in Figs.^Jd-f), where 
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a slightly narrower initial wavepacket was used. As a 
general guide, the peak density and width of the initial 
wavepacket required for the formation of the fundamen- 
tal gap soliton are related as: A"^ ~ \D\/{w^giD) (see 
[sif. eh. 5). While the time scale in Figs. 2Ia,d) is not 
sufficient to see the shape of the wavepacket stabilize, as 
we evolve the wavepacket even further, it is clear that the 
oscillations are gradually being damped as excess atoms 
are removed into the low-density background. 

The spatial extent of the initial wavepacket is required 
to be large 90 lattice sites), so that the momen- 
tum distribution is very narrow compared to the size 
of the Brillouin zone. The acceleration process works 
best with a narrow momentum distribution and hence a 
wide wavepacket, which is broader than the soliton that is 
eventually formed. If the momentum distribution is too 
broad, undesirable effects may occur, e.g. the wavepacket 
can develop a strong asymmetry as it passes through the 
zero dispersion point, as described in Ref. 0. A wide 
initial wavepacket also reduces the maximum density in 
the initial wavepacket used to form a fundamental soli- 
ton since the width of the fundamental soliton's envelope 
near the band e dge is inversely proportional to its ampli- 
tude, A l/w |3ll |. 

For both the ramping and acceleration processes, 
we have checked the momentum components of the 
wavepacket to ensure that there is no significant exci- 
tation of the upper bands. As discussed above, the rate 
at which the lattice can be adiabatically accelerated de- 
pends very strongly on the size of the energy gap. The 
size of the mini-gaps for superlattices with e ^ 1 are 
much smaller than the size of the gaps of a single-periodic 
lattice for the same lattice height. This makes the above 
described techniques for condensate preparation at the 
edges of a superlattice mini-gap rather difficult to use 
due to the large acceleration times. 

Difficulty in accelerating the condensate to the gap 
edge for narrow gaps produced in shallow single-periodic 
optical lattices was also noted in recent experiments 
However, smaller gaps are useful because they poten- 
tially enable greater values of the effective dispersion [see 
Fig. Ola)]. Below, we study an alternative scheme for 
gap-soliton generation in a stationary lattice. Since the 
problems of large time-scales associated with conditions 
of adiabatic ramping and acceleration are avoided in this 
case, the method would be equally useful for low mod- 
ulation height single-periodic optical lattices and super- 
lattices with narrow mini-gaps. 



IV. GAP-SOLITON GENERATION IN A 
STATIONARY LATTICE 

An alternative method for creating a wavepacket with 
the correct quasi-momentum and internal structure to 
produce spatial gap solitons was first suggested theoret- 
ically in the context of nonlinear optics [3. Hol ITU, and 
was recently employed in experiments on weakly coupled 



waveguide arrays and optically-induced photonic lat- 
tices ^2 01- Apphed to matter waves in optical lattices, 
this technique means that, instead of using a moving lat- 
tice to gradually drag a EEC wavepacket to the edge 
of the Brillouin zone, we start with two non-stationary 
wavepackets with opposite momenta (i.e. k — 1 and 
fc = — 1) corresponding to the Bragg reflection condi- 
tion in a stationary lattice. Their interference produces 
a matter wave with the internal structure resembling that 
of the Bloch wave at the gap edge. This is expected since 
the Bloch wave at the edge of a gap is a periodic standing 
wave formed at the Bragg reflection condition and hence 
can be presented as the superposition of identical for- 
ward and backward moving travelling waves. Recently 
this method was used to theoretically demonstrate the 
formation of optical spatial gap solitons in binary waveg- 
uide arrays, the nonlinear optical analogue to our partic- 
ular optical superlattice svste m for a focusing (at- 
tractive) nonlinearity. In Ref. [ij it was found that this 
method allows for the efficient production of spatial gap 
solitons at the top of the mini-gap. Here we apply this 
method to repulsive condensates (and hence the soliton 
formation at the bottom of the mini-gap). 

We assume the initial state of condensate evolution in 
the form of two identical Gaussian wavepackets, V'l (a;, t — 
0) and 4'2{x, t ~ 0), with momenta ki = k and fc2 = —k 
and phases 9i and O2 respectively. These wavepackets 
could be created using Bragg scattering techniques. As- 
suming that both wavepackets are initially centered at 
the origin, their superposition is given by the expression: 

^{x, t = 0) = Ae-^^'/""' [e^'^^+'^i + g-^kx+io^^ ^ 

where A is a constant amplitude and w is a width. By 
setting the correct fc, 61 and 62, the symmetry of ^{x, t = 
0) can be matched to that of the Bloch waves at various 
gap edges. For example, with fc = 1, = ^2 = 0, Eq. Q 
simplifies to ^{x,t = 0) = v4exp(— x^/w^) cos(x), which 
has the spatial structure of the Bloch wave at the bottom 
edge of the first mini-gap of a superlattice potential [see, 
e-g-, Fig. Ela)]. Setting k ^ I, Oi = 37r/2, and 6*2 = 
7r/2 gives ip{x,t = 0) —Aexp{—x'^/w'^)sm{x) which 
mimics the Bloch wave at the top edge of the first mini- 
gap [see, e.g., Fig.l^b)]. 

In what follows we study the generation of matter- wave 
gap solitons in the two cases discussed above. That is, 
(a,) k = 1, 01 = and 62 = i.e. at the bottom of the 
first superlattice mini-gap and (b) k ^ 1, 61 = 3tt/2 and 
62 = 7r/2 i.e. at the top of the first superlattice mini-gap. 
We consider wavepackets containing N = 1000 atoms, 
with the width of the spatial envelope set at w = 30 
which means the initial wavepacket occupies around 30 
wells. With this soliton excitation method, we can use 
wavepackets with small widths that are closer to the fun- 
damental soliton width at the corresponding atom num- 
ber, as a small initial momentum distribution is not cru- 
cial. If the initial wavepacket has a smaller width, then 
we can use a larger nonlinearity, i.e. larger peak density. 
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(a), 



FIG. 5: Interference patterns of two Gaussian wavepackets 
as described by Eq. JSJl with A — 1.3, k = 1, w — 10, and (a) 
01 = 0, 6I2 = 0, and (b) di = 37r/2, 62 = tv/2. Dotted lines 
show the superlattice potential (e = 0.3 and Vo — 1). 



(b) 
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FIG. 6: Linear {gm ~ 0) evolution of wavepackets formed 
by the interference of two Gaussian wavepackets described 
by Eq. in a superlattice potential (e = 0.3 and Vb = 1). 
Number of atoms is A*' = 1000, k — 1, and w = 30. In (a) the 
phases of the two wavepackets are 61 = and 62 = 0. In (b) 
the phases are 9i — 3n/2 and O2 — respectively. 



and still obtain a fundamental soliton |3l| . Using Eq. Q 
as the initial condition, we solve the ID evolution equa- 
tion ^ with the stationary superlattice potential 
by employing a Fourier split-step method implemented 
using the code generator XMDS 13j| 

In the absence of nonlinearity {gio = 0), we expect 
that the initial wavcpacket disperses and spreads out. 
FigureElshows this occurring for both cases (a) and (b) in 
a superlattice with e — 0.3 and Vq — 1. The wavepacket 
at the top of the mini-gap [Fig. EJb)] exhibits slightly 
greater dispersion then the wavepacket at the bottom of 
the mini-gap [Fig.|BJa)], which agrees with the values of 
the effective dispersion coefficient shown in Fig. |31 . 

When a repulsive nonlinearity {giD > 0) is "turned 
on" , the interplay between negative dispersion and pos- 
itive nonlinearity [gi^D < 0) results in the localization 
of the initial wavepacket into a gap soliton at the bot- 
tom of the mini-gap [see Figs.lT^a-c)]. In contrast, at the 
top of the mini-gap the repulsive nonlinearity interacts 
with the positive dispersion {gioD > 0) to accelerate 
the spreading and breakup of the BEC wavepacket [see 
Fig. Efd)]. When an attractive nonlinearity {giD < 0) is 
used instead, the wavepacket forms a gap soliton at the 
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FIG. 7: Evolution of wavepackets described by Eq. with 
repulsive interaction gm = 0.001, A'' — 1000, fc = 1, and 
w = 30. In (a)-(c) 6^1 = and 82 = i.e. the wavepacket 
is at the bottom of the first mini-gap. In (d) Oi = 3n/2 and 
6*2 ~ n/2 i.e. the wavepacket is at the top of the first mini- 
gap. Each of the superlattice potentials have the same period 
and lattice height Vb = 1. The shape of each superlattice 
however varies: (a) e = 0.2, (b) and (d) e — 0.3, (c) e — 0.5. 



(a). 




FIG. 8; Evolution of a wavepacket formed by the interference 
of two Gaussian wavepackets as described by Eq. Q in a 
superlattice potential with Vb = 1 and e — 0.3 for A'^ = 1000, 
k = 1, and w = 30. In contrast to Figs. |7Ib) and (d), we use 
an attractive nonlinearity, gm ~ —0.001. In (a) 9i = and 
6*2 = i.e. the wavepacket is at the bottom the first mini-gap. 
In (b) 61 = 3n/2 and 02 = 7i"/2 so that the wavepacket is at 
the top of the first mini-gap. 



top of the mini-gap while the one at the bottom breaks 
up more quickly [see Fig. In the three cases shown 
in Figs. [Tfa-c), the fraction of atoms radiated from the 
wavepacket during the process of gap-soliton formation 
is less than 10%. In addition, the localization occurs 
over time intervals that are almost an order of magni- 
tude shorter than those for the transient dynamics of the 
wavepackets in the moving lattices. 

In general, the interference method of gap-soliton gen- 
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FIG. 9: Left: Density profiles of the pulses shown in Fig- 
uresUJa)-(c). The initial conditions are identical in each case, 
but the superlattice potential varies with (a) e = 0.2, (b) 
£ = 0.3 and (c) £ = 0.5. In each case the lattice height is 
Vb = 1. Right: Exact stationary gap soliton solutions ob- 
tained by numerically solving the time-independent version 
of the full model equation (1) for the same superlattice pa- 
rameters as on the left. 



eration allows for more precise control of the initial 
wavepacket widths, uninhibited by any requirements on 
the initial momentum distribution. This allows us to per- 
form an accurate analysis of the effects of dispersion on 
the BEC localization process in the superlattice poten- 
tial by creating an initial wavepacket at the band-gap 
edge with exactly the same parameters for each super- 
lattice. The effects of the dispersion control are demon- 
strated in Figs. |3a-c), which show the gap-soliton gen- 
eration in a mini-gap for three different shapes of the 
superlattice potentials and the same repulsive nonlinear- 
ity {giD = 0.001). In each case, the lattice period and 
height are kept constant, as are the parameters of the 
initial wavepacket. The superlattice parameter e varies 
and with it the relative well-depths of the small and large 
wells in the superlattice potential. With growing s, the 
generated gap soliton becomes narrower and the aver- 
age peak density increases. The dramatic change in the 



shape of the output soliton with variations in e is due to 
the variations in the effective dispersion with e, as seen 
in Fig.Etb). 

Figure O shows that as e ^ 1 the magnitude of the 
effective dispersion at the edges of the superlattice band- 
gap decreases due to the increasing size of the mini-gap. 
Consequently, the peak particle density required to form 
a soliton with the same number of atoms needs to increase 
to compensate. In Fig.|51we show the density profiles of 
the pulses shown in Figs. I3a-c) compared with the sta- 
tionary fundamental soliton solutions obtained numeri- 
cally by solving the time-independent GP equation for 
identical superlattice parameters. Because of the tran- 
sient oscillations in the peak density of the dynamically 
generated gap solitons, we use the average value over 
time for this comparison. In each case, the matter-wave 
chemical potential is around the value fi = 1.0, which is 
close to the bottom edge of the mini-gap for each of the 
three values of e. 



V. CONCLUSIONS 

In this paper, we have investigated an effective method 
of generating matter-wave gap solitons - nonlinear local- 
ized states of a repulsive Bose-Einstein condensate in a 
stationary optical lattice. We demonstrated that, by in- 
terfering two matter waves, one can create a wavepacket 
located at the gap edge with the appropriate internal 
structure suitable for gap soliton formation. The ini- 
tial wavepacket prepared in this fashion can be better 
matched to the shape of the fundamental gap soliton, 
which leads to a more efficient and faster soliton forma- 
tion compared to the standard technique of lattice ac- 
celeration. This method would be especially useful for 
small energy gaps, e.g. in the case of superlattice mini- 
gaps and for extremely shallow single-periodic optical lat- 
tices, where preparation of the condensate wavepackets 
at a gap edge by using lattice acceleration is experimen- 
tally difficult and time consuming. Using the interference 
method to simulate the gap soliton formation in double- 
periodic optical superlattices, we have shown that it is 
possible to fine-tune the effective dispersion of matter 
waves by changing the shape of the superlattice rather 
than its height or periodicity, and hence tailor the prop- 
erties of the BEC solitons. 
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